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r~| . We study in detail the effect of the finite current quark mass on chiral symmetry 

Qh| breaking, in the framework of truncated Coulomb gauge QCD with a linear con- 

{^' fining quark-antiquark potential. In the chiral limit of massless current quarks, the 

^ . breaking of chiral symmetry is spontaneous. But for a finite current quark mass, 

some dynamical symmetry breaking continues to add to the explicit breaking caused 
^^.^ by the quark mass. Moreover, using as order parameter the mass gap, i. e. the quark 

^ . mass at vanishing moment or the quark condensate, a finite quark mass transforms 

'sj" I the chiral symmetry breaking from a phase transition into a crossover. For the study 

of the QCD phase diagram it thus is relevant to determine how the current quark 
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^\j ■ mass affects chiral symmetry breaking. Since the current quark masses of the six 

standard flavours u, d, s, c, b, t span over five orders of magnitude from 1.5 MeV 

^^ ' to 171 GeV, we develop an accurate numerical method to study the running quark 

^P ■ mass gap and the quark vacuum energy density from very small to very large current 

quark masses. 
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1 Introduction 



Chiral symmetry breaking has been studied in detail with chiral invariant and 
confining quark models in the chiral limit of a vanishing current quark mass 
niQ. For a finite current quark mass mo ^ 0, studies exist with an approximate 
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confinement [1,2,3,4] or with a quadratic confining potential [7,6,5], but very 
few studies have been performed [8,9] with an exactly linear confining poten- 
tial. Since the current quark masses of the six standard flavours u, d, s, c, b, t 
span over five orders of magnitude from 1.5 MeV to 171 GeV, here we develop 
a new numerical method to study the quark mass gap and the quark vacuum 
energy density, with a linear exactly confining potential, from very small to 
very large quark masses. 

Notice that even in the chiral limit of rriQ = 0, the quark has a constituent 
running mass m{p) function of the momentum p, solution of the mass gap 
equation (equivalent to the Schwinger-Dyson equation) for the quark. Recently 
we have shown how to measure in the excited hadron spectra the running mass 
m{p) [10]. The chiral invariant and confining quark models have also been 
applied to phase studies at finite temperature T and chemical potential /i 
[11,12,13,14,15,16,17,18]. A finite quark mass is relevant both for the study of 
hadrons which have been investigated for decades, and for the study the QCD 
phase diagram which will be explored in the future at RHIC, LHC and FAIR. 
In the phase diagram, a finite current quark mass niQ affects the position of 
the critical point between the crossover at low chemical potential fi and the 
phase transition at higher fi. Moreover the current quark mass affects the QCD 
vacuum energy density S, relevant for the dark energy of cosmology. This all 
occurs in the dynamical generation of the quark mass m{p). While the quark 
condensate {ipip) is a frequently used order parameter for chiral symmetry 
breaking, the mass gap, i. e. the quark mass at vanishing momentum m(0) 
is another possible order parameter for chiral symmetry breaking. However, 
due to technical difficulties, m(0) has not been computed in detail previously 
in confining and chiral invariant quark models. 

Here we study in detail how a finite current quark mass ttiq ^ affects the 
dynamically generated quark mass m{p) . We utilize the linear confining poten- 
tial for the quark-antiquark interaction, in the chiral quark model or Coulomb 
Gauge quark model, including both confinement and chiral symmetry. While 
this model, in the framework Coulomb gauge Hamiltonian formalism is not 
yet full QCD, it is presently the only model able to include both the quark- 
antiquark confining potential and quark-antiquark vacuum condensation. Im- 
portantly, since our model is well defined and solvable, it can be used as a 
simpler model than QCD, and yet qualitatively correct, to address different 
aspects of hadronic physics. Is is adequate to study the QCD phase diagram 
microscopically [11,12,13,14,15,16,17,18]. We scan the current quark masses, 
from the light quarks to the heavy quarks, computing the running quark mass 
mij)) with detail, including the infrared limit of ?7i(0), i. e. the mass gap. 

In Section II we derive in detail the mass gap equation. In Section HI we review 
the numerical difficulties of this non-linear integral equation, with cancelling 
infrared divergences. We solve the mass gap equation with a new numerical 



method, dedicated to determine in detail the difficult infrared region of small 
momentum p ~ 0. In Section IV we discuss our results and conclude. 



2 Our framework 



2.1 Possible link to QCD. 



Our framework can be approximately derived from QCD, in two different 
gauges. In Coulomb gauge V • A(x, t) = the interaction potential, has been 
derived by Lee, [19], and by Szczepaniak and Swanson [20,21]. In the present 
study we address the quark fields only and thus Vj reduces to the quark part 
of the density- density term, 

Vi=+-g^ /"dxrfyJ^-V^(x)TXx)(x,a| x 

(V . D)-i(-V2)(V • D)-V, 6):r^t(y)TV(y) 

(1) 

The covariant derivative in the adjoint representation D = V — gA, and 
J' = Det[V-D] contribute to the density-density interaction, which is expected 
to be confining in QCD. 

Another approximate path from QCD to our model considers the modified 
coordinate gauge of Balitsky [22], A(0,t) = , x ■ A(x, t) = and in the 
interaction potential for the quark sector. 



Vt= / d-^x 



ip^'ix) (mo/3 -id-V) tpix) + -g"^ / d^y 



2^ 

i^{x)r'-irH^){A''^{x)Al{y)) i^{y)Y^ijiy) + ■■■ (2) 



retains the first cumulant order, of two gluons [23,24,25] (^^(^^(x)^^^/)) and 
this also results in a simple density-density harmonic effective confining inter- 
action. As in QCD, this only has one scale. 

Thus our framework is similar to an expansion of the QCD interaction, trun- 
cated to the leading density-density term. Moreover, to address phenomenol- 
ogy where the meson spectrum fits in linear Regge trajectories, one also needs 
to assume that the confining quark-antiquark potential is a linear potential. 
Notice that the short range Coulomb potential could also be included in the 
interaction, but here we ignore it since it only affects the quark mass through 



ultraviolet renormalization [26], which is assumed to be already included in 
the current quark mass. While this is not exactly equivalent to QCD, our 
framework maintains three interesting aspects of non-pertubative QCD, a 
chiral invariant quark-antiquark interaction, [27,28,29,30,31,32] the cancel- 
lation of infrared divergences [7,6,5] and a quark-antiquark linear potential 
[8,9,20,33,34,35,36]. 



2.2 Deriving the mass gap equation 

We derive the mass gap equation, where constituent quarks acquire the con- 
stituent mass m{k) [37] in the true and stable vacuum, solving the Schwinger- 
Dyson equation for the quark propagator, 

p—k 

S^\p)=So^\p) - O. (3) 



We utilize the truncated Schwinger-Dyson equation at the Rainbow level, 
where the dotted line represents the same density- density interaction Vj re- 
sulting identically from the truncation of Coulomb gauge QCD or of Balitsky 
gauge QCD. This leads to the same mass gap equation and quark dispersion 
relation obtained assuming a quark-antiquark ^Pq condensed vacuum, com- 
puting the vacuum energy density with the Hamiltonian, and minimizing the 
energy density. The relativistic invariant Dirac-Feynman [31] propagators S, 
can be decomposed in the quark and antiquark Bethe-Goldstone propagators 
[37], 



S{ko,k) 



/c — m + ie 



(4) 



ko — vW+m? + ie 



Y^UsulP 



— /cq — a/P + jt? + ie 



j:^svif3 



where m is the quark mass and where the quark spinor Ug and antiquark spinor 
Wo are, 




k-a-f5 



k- a'j5 



UsiO) 



(5) 



where S = m/y/W+lv? is a function of the quark mass. 



Importantly, in the free propagator, the correct quark propagator in the non 
condensed vacuum, the quark mass m is equal to the current quark mass niQ. 
And it is this current quark mass tuq which effects in the current quark mass 
we study in great detail. However when chiral symmetry breaking occurs, m is 
not determined from the onset. In the physical vacuum, the constituent quark 
mass m{k), is a function of the momentum, a dynamical solution of the mass 
gap equation. 

Replacing the propagator of eq. (5) in the Schwinger-Dyson equation and 
projecting it with the spinors, we get the mass gap equation and the quark 
dispersion relation. 



kk-cx + moP-j ^^s^yik - k')Siko\ k') \ v,.{k) (6) 

{r dk ' r/'^k' ~ 1 

kk-cx + mo/3-J -^T^^Vik - k')S{k^\ k') \ Us{k), 

where the usual notation for Dirac matrices is assumed. Writing the running 
mass in terms of a sine and a cosine of </?(A;) = arctan ^^, the chiral angle, 



S{k) = sm^{k)- ""^^^ 



Jk'^ + m{ky 

C{k) = cos ip{k) = ^ , (7) 

Jk"^ + m{ky 

the mass gap equation and the quark energy are. 



= +S{p)B{p)-C{p)A{p) (8) 

E{p) = +S{p)A{p) + C{p)B{p) (9) 

where the propagator functions A and B, respectively replacing the quark 
mass m and quark momentum |p| in the one-loop dressed quark propagator 
of eq. (4) are. 



1 r d^h ~ 
B{p) =P+^J ^^^(P - k)(P ■ k)C{k) . (10) 



Equivalently to solve the non-linear integral mass gap equation (8), we can 
alternatively minimize the vacuum energy density per unit volume, 

^ = -f / -0J3S{p) [A{p) + mo] + C{p) [B{p) + p] (11) 



where g = Nf Ng Nc is the degeneracy factor counting the number of different 
but degenerate quarks. A^^ = 2 is the number of spins and A^c = 3 is the 
number of colours. A'^^ is the number of degenerate flavours, but since each 
quark has a different current quark mass uiq one should compute separately 
the vacuum energy difference for each quark flavour. 



2.3 The mass gap equation for a linear confining potential 



Notice that in the case of a linear potential, divergent in the infrared limit of 
large r, the Fourier transform needs an infrared regulator //. A possible form 
of the linear potential, 

V(r) = -a- ~ -- + (Tr-— r2 + --- (12) 

jj, 11 2 

corresponds, in the limit of small infrared regulator //, to a model of linear 
confinement where the quark also has an infinite binding energy — -. While 
other infrared regulations can be used for the linear potential, say V{r) = 
ar e~^ *" which has no infrared divergent binding energy, the infrared divergent 
constant of Eq. (12) is exactly cancelled in the mass gap equation by the factor 
in the integrand [S{k)C{p) — S{p)C{k)k ■ p\. The potential in Eq. (12) has a 
simple three-dimensional Fourier transform. 



and this is the momentum space potential frequently utilized to account for 
linear confinement. 

The integrals in the angular variable u of Eq. (10) are. 
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{2kp)n [{k - py + i2^][{k + p^ + fi^] 2 
We find for tlie propagator functions A and B, 
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(15) 



leading to the mass gap equation in the two equivalent forms of a non-linear 
integral functional equation, 



a r dk 



0=pS{p) -moC{p) 

p^ J In 
^ 

IA{k,p,^l)S{k)C{p)-IB{k,p,^l)S{p)C{k)] 

and of a minimum equation of the energy density, 



(16) 



2n J 27r 



dk 



2p^C{p) + 2p^moS{p) + a X 



—lA{k,p,^i)S{k)S{p) + lB{k,p,^x)C{p)C{k) 

2tx 



(17) 



Eq. (16) can be rewritten as a fixed point equation for the quark mass function 
m(/c). 



a f dk I^ikyP, fi) m(k)p — lB(k,p, fi) m(p)k 
P^ I 2vr Jk^ + m{ky 



(18) 



Since the potential has an infinite constant independent of the mass, in the 
variational equation we search for the extremum of the energy difference 



8-£, = ^ fM2p'[C{p) - Cob)] + 2p^mo[S(p) - So{p)] (19) 

ZTX J ZTT L 


—I^{k,p,fi) [S{k)S{p) - So{k)So{p)] 
Zn 



+lB{k,p,fi) [C{p)C{k)-Co{k)Coip)] 



where Sq is constant, and where we use the subindex q when the mass m^p) 
is substituted by the constant current mass niQ. These two Eqs. (18) and (19) 
constitute the main object of our study. 



3 Solving the mass gap equation variationally with rational ansatze 

3. 1 Accurate numerical cancellation of infrared and ultraviolet divergences 



We use both Eq. (18) and the minimization of Eq. (19) to find the dynamical 
quark mass m{k), but first we must regulate correctly their divergences. The 
infrared divergences are present in the term pk/[{p — k)"^ + /i^], infrared diver- 
gent in the limit of a vanishing regulator /i — )■ 0, which is present bot in the 
functions /^(/c,p, /i) and lB{k,p,fi). We must show that this infrared diver- 
gence is cancelled both in the fixed point Eq. (18) and in the energy density 
Eq. (19). In what concerns the fixed point Eq. (18), while the denominator 
diverges quadratically in l/{k —pY, the numerator m{k)p — m{p)k is of order 
(k — p) and thus the integrand diverges like l/{p—k) only, and it's integral has 
finite principal value. In the energy density difference the infrared divergence 
also cancels, since the numerator common to the infrared divergent terms. 



Sip)S{k) - So{p)So{k) + C{p)C{k) - C,{p)C,{k) 
= -\^y\p)~v',\p)]{k-pf (20) 

-\\v\pW\p) - ^oiP)^M]{k - Pf + o{p - kY 

is then of order {k — p)"^. 

In the ultraviolet part of the integrals, while each sub-term in the propagator 
function integrands lA{k,p,fj,) and lB{k,p,fj,) is divergent, the actual sum is 
ultraviolet convergent, since in the limit of large k we have. 



oo „2n+l 



n=l 

oo 




fcjB(fc,P,/x) = E4n|;::^)^, (21) 

n=l 

and thus the integrals in k have ultraviolet integrable integrands decaying 
like p^/k^. Also, the ultraviolet divergence in the kinetic terms of the energy 
density 2p^C{p) + 2p'^moS{p) cancels due to the difference with 2p^Co{p) + 
2p'^'mQSo{p) if the mass difference mlk) — mo vanishes sufficiently fast in the 
ultraviolet. 

To address correctly both the infrared and ultraviolet divergences of the in- 
tegrals in k of Eqs. (18) and (19), we divide the integral in two sections, the 
infrared one ioi < k < 2p and the ultraviolet one for 2p < k < oo. In 
the infrared region we compute the respective principal value, performing a 
symmetric sum centred in k = p, maintaining a very small regulator fi just to 
cancel automatically the contribution oi k = p. In the ultraviolet region we 
use the change of variable [33] of Adler and Davis k — > x/{l — x) with Jacobian 
1/(1 — x^ and with integration between x = 2p/{l + 2p) and 1. The change 
of variable in the ultraviolet transforms, say an integral of rational functions 
1/(1 + k)^ into the integral of polynomials (1 — x)"~^ and thus it is adequate 
to the integral of rational functions as we have here. Our numerical integrals 
in /c of a generic integrand X singular in p are thus computed in the form, 

'" rl (t^) . (22) 
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where each of the two numerical integrals can either be computed with a 
rectangular, trapezoidal or Simpson sum or with the gauss quadrature method. 

I one would discretize the quark mass m^p) in a series of momenta Pi, then 
the integrals of Eq. (22) loses accuracy. Finite differences would require many 
interpolations, both in the infrared end and in the ultraviolet end, since the 
principal value requires requires that k has many summation points smaller 
than p and many other larger than p. Moreover the correct integration of the 
integrand in the ultraviolet large k limit, where the integrand behaves like 
f I j , also requires an integration extending beyond any value of p. To solve 
this problem, we utilize a well defined ansatz for m{k), formally describe the 
parametrized as. 



m{p) = m{p; ci, C2, ■ ■ ■ , c„) . (23) 

and this allows the numerical summation for the integrals in any point of the 



integration domain. 

In what concerns the numerical convergence to the solution for m{k), the 
fixed point equation is relatively unstable, particularly in the infrared region 
of p ^ 0, when we are searching for the vacuum groundstate. Notice that the 
mass gap equation had not only one, but an infinite tower of solutions [38]. The 
excited solutions dominate the fixed point iteration, converging to the larger 
eigenvalue of the iterated matrix, because they minimize the denominator 
\/k'^ + m?. Previously in the literature, the fixed point method was provided 
with extra stability with two different methods, Adler and Davis used a cubic 
equation and relaxation, to select the best solution [33]. Bicudo and Nefediev 
quasi-linearized the fixed point equation and selected the desired eigenvalue, 
corresponding either to the stable vacuum or to excited, false vacua [38]. Thus 
they were able to find both the stable vacuum and the excited false vacua. 
But these works have not yet determined in detail m(0), since this demands a 
very large numerical precision, and since most previous authors have focused 
in computing the function S{p) which in the infrared region is 5(0) = 1 
regardless of the actual finite value of m(0). Importantly, the present technique 
of minimizing the energy density directly tends to the right solution, which is 
the groundstate vacuum. 

Interestingly, the variation of the ansatz parameters ci, C2, ■ ■ ■ , c^ of the energy 
density of the vacuum £ = S{ci,C2,- ■ ■ ,Cn) utilized in minimization codes 
with gradient method, utilizes the fixed point equation. Computing the partial 
derivatives of the energy density we get. 



dS(ci) /■ , S£ d(p 
dp 



dci J 5(p dci 

= -f J 'A^2p')n(p; c.) . ^% Z'"^^^-' . (24) 

27r J 2% p^ + m^{p] q) oci 

where 7l{p; q) is the right hand side of the mass gap Eq. (8), 

TZ{p; Ci) = +S{p; Ci)B{p; q) - C{p; Ci)A{p; d) , (25) 

utilized in the fixed point Eq. (18). 

3.2 Choosing variational ansatze for the quark mass m{p) 



To guide our choice of ansatze m(p; Ci, C2, • • • , c„), we first notice that the series 
expansion for I a and Ib in Eq. (21), also apply when k -H- p. When replaced 
in the integral of the fixed point equation (18) , the series suggest that a series 
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expansion of ?7i(p) should only have even terms, i. e. m{p) should be a function 
of p^. m{p) should also be a finite function since our integrals are finite. 

In what concerns the asymptotic ultraviolet tail of the integral of the fixed 
point equation (18), there are two different limits we can address. In the case 
of a large current quark mass mo, k'^T^ 2 interpolates between 1 in the infrared 
region of the integral and ^ in the ultraviolet region of the integral. Using 
these approximation, the components of the integral are analytical. 



,3 



-paA;,p,;z) = -, (26) 
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—klBik,p,fi) = —- , 
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dk 
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klBik,p,n) = - 
'^ 
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//i TV 

dk ^ ,, , -PAt + (P^ + /^^) arctan ^ 



00 



J ^^^'^^'^'^^ 



TT/X 



and thus adding the respective components we find that in the infared 
dominated approximation the integrals cancel, while in the ultraviolet domi- 
nated approximation the integral with ^ produces an ultraviolet behaviour 
of m{p) — mo — > ^^. Thus in the case of large ttiq we expect that m{p) — niQ 
decays in the ultraviolet like 1/p^. 

In the case where mo ~ 0, assuming then that in the large p limit the dy- 
namical quark mass vanishes sufficiently fast, the fixed point equation leads 

to, 



%2 



, , a f dk 1 ,^ /,N 

^ ' P' i 27r ^k^^m{kf P' ^ ' 



Aa r dk k'^m{k) 
~^ P^i 27r Jk^ + m{ky 



and, providing m{p) decays faster than l/p"^ for a finite integral, this decays 
in the ultraviolet like 1/p^. 
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The simplest possible ansatz for m^p) — rriQ we may have, function of p"^, 
and encompassing both the behaviour in 1/p^ for a large current quark mass 
mo and the behaviour in 1/p^ for a small current quark mass is the rational 
function, 



Asip) = ^ \^ J- (28) 

Co + C2P^ + C4P* 

This ansatz is a Pade approximant, and to check whether our simple ansatz is 
sufficient, it is convenient to check that the next Pade approximant, a more 
flexible rational ansatz with two more parameters. 



AJv) = ^ + ""'P' (29) 

leads to the same result. In both ansatze of Eq. (28) and Eq. (29) we assume 
that the parameters are positive. While Eq. (28) is a decreasing function, 
Eq. (29) may have a different behaviour at the origin, either with an initial 
increase, or with a steeper decrease, an thus it has room in it's parameter 
space to verify if the ansatz of Eq. (28) is close to the correct solution of the 
mass gap equation. 

We also check numerically that ansatze with steeper infrared behaviours, in- 
cluding in the denominator terms like a Ci/c or a C-i/k would not improve 
the solution since the best solution would have ci = c_i = 0. A better ansatz 
than AsIp) is A5{p), however the improvement of the solution is very small, 
almost invisible to the naked eye in graphics. The partial redundancy between 
the numerator and denominator parameters of A^lp) already slows the con- 
vergence to the minimum of the energy density S. Thus an ansatz with more 
parameters than A^lp) is not necessary. The only ansatze we adopt here are 
the ones of the rational functions A3{p) and A^lp). 



3.3 One loop results for the large current mass m^ » -y/a limit 



We now compute the first iteration of the fixed point method starting with 
m{k) = rriQ. In the simple case case of a constant mass m{k) = mo, we 
can compute the integral in the fixed point equation with a large precision. 
Defining the mass difference V, 

T>{p) = m{p) — niQ , (30) 
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we compute 'P(p) in the case a constant mass rriQ is used in the integrand, 
mocr JdklA{k,p,fi)p-lBik,p,iJ,)k 



We get an accurate result for the integral Vq{p), with a numerical integration 
decomposed according to Eq. (22). 

This provides a good quark mass solution to m^p) = rriQ + V{p) whenever the 
current quark mass is large, i. e. when mo » T>{p). In that case the integral 
Vq^p) only needs to be computed once, since this one- loop approximation is 
already excellent. 

Moreover we can rescale in mo, and then with a single computation we get 
m{p) for for any constant mass mo- Denoting k = k/rriQ and so forth we get. 



a f dklA{k,p,^)p- lB{k,p,fi)k 

^o(p) = zr^ 



niQ mo 

OO 

1 f dklAik,p,fi)p-lB{k,p,^)k 

^f"* = ? / 2^ 7PT7 ■ '''* 

Thus we only need to compute the dimensionless function J^{p) and this will 
produce the dynamical quark mass for any current quark mass m,Q larger than 
the typical scale ^/o^ of our problem, 

(J p 

mip) ^ mo H J^i — ) . (33) 

mo mo 



The solution of the integration, obtained with the simplest numerical rect- 
angular sum in both integrals, but needing 10^ integration points at least, is 
represented in Fig. 1. 

To be able to perform an accurate integration both in the infrared and the 
ultraviolet, it is necessary to have an analytical function. We obtain an analyt- 
ical function by fitting the function J^{p) , and we get an excellent fit already 
with the very simple ansatz Az{p) of Eq. (28) for the parameter set, 

Co = 4.7664, C2 = 3.4762, a = 0.0000, (34) 

with almost no graphically visible difference in Fig. 1 from the ansatz A^i^p) 
with two more parameters of Eq. (29) for the parameter set. 
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Fig. 1. Dimensionless mass difference function J-{p) computed in the limit when 
the current quark mass rriQ is large. The dots show the numerical integral of Eq. 
(32) and the two almost overlapping solid lines show our fits with the two different 
rational ansatze of Eq. (28) and Eq. (29). 

do = 4.7250, ^2 = 4.9818, d^ = 0.8392, 4 = 0.0000, na = 0.2626, (35) 

and both fits confirm the 1/p^ decay of the mass in the ultraviolet, for large 
current quark masses. 

Now we can include the large quark limit scaling in rriQ, and also apply to the 
quark mass difference T>q{p) the ansatz ^3(p) of Eq. (28), and in this case the 
ansatz parameters should be respectively. 



Co = 4.7664 mo , C2 = 3.4762 m(^^ , C4 = 0.0000 m(^^ . 



(36) 



We further compute the vacuum energy density difference for a large current 
quark mass niQ. Although the mass difference T>q{p) in Eq. (32) decreases like 
??7,o~^ in the large current quark mass mo limit, a dimensional analysis oi S — Sq 
shows that possibly it does not vanish. Since the energy difference is finite and 
independent of the infrared cutoff //, then our only dimensionfuU parameters 
are the string tension a (scaling like a mass square) and the current quark 
mass rriQ (scaling like a mass). Because the vacuum energy density per volume 
scales like a mass to the fourth power, a large mass expansion may have terms 
of the form rriQ'^, rriQ^a, o"^, a^mQ~^, ■ ■ ■ However the first term in this series 
clearly vanishes in the vacuum energy difference £ — £q. Then the question is 
what is the first non-vanishing term in this expansion. 
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mo 


IqAE-So 


Co 


C2 


C4 


0.000100000 


-0.473813 


6.24900 


26.7910 


17.5059 


0.000316228 


-0.478133 


6.24287 


26.7090 


17.3392 


0.00100000 


-0.491581 


6.22441 


26.4168 


16.8674 


0.00316228 


-0.532245 


6.17597 


25.4076 


15.8086 


0.0100000 


-0.649468 


6.01623 


23.2517 


12.0965 


0.0316228 


-0.958565 


5.68682 


18.5897 


6.38068 


0.100000 


-1.62048 


5.22374 


12.5110 


1.50433 


0.316228 


-2.58117 


5.04452 


7.22212 


0.051756 


1.00000 


-3.27613 


7.03428 


3.06500 


0.000000 


3.16228 


-3.46213 


17.1365 


1.04266 


0.000000 


10.0000 


-3.47122 


51.7243 


0.336172 


0.000000 


31.6228 


-3.39336 


153.536 


0.113085 


0.000000 



Table 1 

The results for the ansatz parameters of A3{p) for m{p) — mo obtained with our 

minimization code. Ah results are in dimensionless units of cj = 0.19 GeV^ = 1. 



To answer this question, we expand the vacuum energy density in a ^^ series, 
similar to the variation in Eq. (24), but now based in the expansion of the 
quark dynamical mass. 



m{p) 



mo H nil — 

nio Knio 



(y f P 

— 3^^ — 
TTio Vrrao 



+ ■ 



(37) 



and then we get 



9 [ dp \ o„.2 



27r y 27r ^ ^ 



' n ( \ ^ Q I \ ^ pSCo{p)+mo6So{p) , , . . 

P^o{P) + moSoip) H {m{p) - mo J 



+ 2 



lp6^C{p) + mo6^So{p) 



5mo^ 



6nio 



{m{p) - mo) + 



(3J 







27r 



So{p)So{k) + 2^-^So{k) {m{p) - mo) 
omo 



+lB{k,p,n) 



Coip)Co{k) + 2^-^^Coik) im{p) - mo) + 
omo 



The zeroth order term vanishes when we perform the difference of the vacuum 
energy densities S — Sq. Then the first order term in the kinetic energy density 
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mo 


104^ 


do 


d2 


d^ 


4 


^2 


0.000100000 


-0.473908 


6.170 


31.9 


22. 


14. 


0.5 


0.000316228 


-0.478221 


6.163 


32.3 


23. 


15. 


0.6 


0.00100000 


-0.491647 


6.147 


32.2 


24. 


16. 


0.6 


0.00316228 


-0.532266 


6.09 


33.1 


30. 


19. 


0.9 


0.0100000 


-0.649397 


5.9705 


29.582 


28.37 


12.83 


0.871 


0.0316228 


-0.958497 


5.7202 


25.358 


32.40 


7.4897 


1.301 


0.100000 


-1.624761 


5.369 


11.875 


5.05 


0.0000 


0.149 


0.316228 


-2.59104 


5.467 


6.720 


2.08 


0.0000 


0.205 


1.00000 


-3.27722 


7.5349 


4.4398 


0.8105 


0.0000 


0.2546 


3.16228 


-3.46240 


16.665 


1.2612 


0.0105 


0.0000 


0.0104 


10.0000 


-3.47148 


49. 


0.6 


0.00 


0.0000 


0.00 


31.6228 


-3.76155 


92. 


0.2 


0.0000 


0.0000 


0.0000 



Table 2 

The results for the ansatz parameters of A5{p) for m{p) — 
minimization code. All results are in dimensionless units of 



niQ obtained with our 
a = 0.19 GeV2 = 1. 



also vanishes since the kinetic energy density is minimized by m^p) = niQ. 
Thus the leading term is of second order in the kinetic energy density, and of 
first order in the potential energy density, and both are proportional to o"^. 
Using the intermediate steps. 



pS^Co{p)+mo6^So{p) 



p" 



a 



dk 

2^ 



■ 26Soip) 
. SniQ 



{p^ + niQ^) 



lA{k,p,fi)So{k) + lB{k,p,fi)Co{k) 



25Coip) 
drriQ 



p-* 



a 



D„ (^-\ 



{p^ + mo^f^^rrio \moJ 
and changing variable to the dimensionless p, we get, 



(39) 



£ — So = a 



2-g [ dpp\J^{j>)f ( a- 



2tx J 27r(p2 + i): 



mo" 



(40) 



Finally using the ansatz Az{p) of Eq. (28), with the parameter set of Eq. (34) 
this results in 
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Fig. 2. Plots of our numerical solution with the ansatz A3{p), as a function of the 
current quark mass mo: (a) the vacuum energy density shift S — Sq, (b) parameter 
Co , (c) parameter C2 , (d) Parameter C4. The dots show our numerical solution, the 
solid line is the large itiq limit obtained with Eqs. (36) and (41), and the vertical 
dot-dashed lines represent the current masses of the quarks u, d, s, c, b, t. 



S-Sn 



3.47701 X 10~^ a^g . 



(41) 



4 Results 



Utilizing our ansatze of Eq. (28) and Eq. (29), we may now compute with 
great accuracy the integrals of Eq. (19) which now are a function of the ansatz 
parameters, and apply a standard minimizing code to determine the optimal 
parameters. We use 1000 x 1000 integration points and up to 40 decimal 
digits in order to be able to find a convergence of the method in the case of 
the ansatz A^lp), due to the partial redundancy of the parameters. Then we 
also minimize the energy starting from different randomly generated initial 
values for the parameters. In Tables 1 and 2 we only show the digits which are 
stable, in the sense that they do not depend on the initial values. Notice that 
the energy density obtained with the two different ansatze differ only by a 
few per mil and that the ansatz of Eq. (29) already exhibits some redundance 
in the parameters. This shows that the ansatz of Eq. (28) is already quite 
accurate for the parametrization of the quark mass m{p). 
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Fig. 3. Plots of the mass gap as a function of the current quark mass mo: (a) the 
mass gap difference m{0) — mQ, measuring the amount of generated dynamical mass, 
(b) the mass gap m(0) . The dots show our numerical solution, the solid line is the 
leading order obtained when itiq — t- oo, and the vertical dot-dashed lines represent 
the current masses of the quarks u, d, s, c, b, t. Notice that the dynamical mass 
generation has a maximum for finite quark masses close to the strange quark mass. 
All results are in dimensionless units of cr = 0.19 GeV^ = 1. 

Unlike the fixed point method, converging quite fast (with a single iteration) 
for large current quark masses rriQ » -\/o", the variational method converges 
faster for small current quark masses. Thus both methods are complementary. 
In Tables 1 and 2 and in Fig. 2 we show the results of our minimization for 
two ansatze and for different current quark masses spanning over five orders of 
magnitude. In Figs. 3 and 4 we illustrate the mass difference at the momentum 
origin ?Ti(0) — rrio, the mass gap 771(0) and the running mass m{p) for different 
current masses mo- 

Now that we have an excellent and simple ansatz for the running quark mass, 
we may compute the regularized quark condensate and the quark dispersion 
relation. The quark condensate (ipip) is another possible order parameter, to be 
compared with the other order parameter we computed, i. e. the quark mass 
at vanishing momentum m(0). The quark condensate is computed from the 
one-loop quark propagator functions in Eq. (10), and it is ultraviolet divergent 
for finite quark masses. Thus we regularize the quark condensate, subtracting 
the quark condensate for a current quark. 



(^7/>) - (^^)o = -^ I k'dk 



m{k) 



rriQ 



^k^ + m{ky Vk'' + mo\ 



(42) 



The one quark dispersion relation E(p), defined in Eq. (9), is relevant for the 
boundstate equation of mesons or of baryons. For instance, in the instanta- 
neous Salpeter equation, a hamiltonian H = Eg + Eg + Vgg can be defined 
for mesons (actually the hamiltonian is a matrix [37] including negative and 



D(p) 
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P 



Fig. 4. The quark mass function difference I'(p) = m{p) — mo measuring 
the extent of dynamical mass generation, is represented with increasing num- 
ber of dashes per curve for five different current quark masses niQ with values 
10"^, 10-^ 10-1/2^ 10°, 10^/2^ in dimensionless units of a = 0.19 GeV^ = 1. 

positive energy components). The dispersion relation E{p) is infrared diver- 
gent due to the infinite constant of quark- ant iquark potential detailed in Eq. 
(12). In momentum space, this leads to an infinite Dirac delta in the integral 
present in E{p). We can regularize the numerical integral, subtracting a term 
to cancel the integrand when k = p, a. term that we add back analytically, 



°° 77 

E{p)=pC{p)+moS{p) + - \l^{k,p,f,) [S{k)-S{p)]S{p) (43) 



+lB(k,p,fi) [C{p)-C{p)]C{k)+ lA{k,p,^i)S\p) + lB{k,p,fi)C\p) 



in particular the integral of /a 5"^ -|- IbC^ is analytical thanks to Eq. (26), 
and in the limit of a vanishing infrared regulator /x — )■ this term, that we 
subtracted and must now add back to the quark energy, reduces to a simple 
infrared divergence plus a finite term. 



^ J ^lA{k,p,f,)S'{p) + lB{k,p,f,)C\p) ^ 



a 2a 



p 



2/i TT p2 — ■m(py 



(44) 



Importantly, this divergence is physically irrelevant since, in the hamiltonian 
of any colour singlet hadron, the sum of the divergences of the quark and of the 
antiquark energies cancel with the infrared divergence of the quark-antiquark 
potential detailed in Eq. (12). Finally for the purpose of future computations 
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mo 


^ ~{^(>1p) + {■4>■4))o\ 4 


eo 


ei 


63 


es 


0.000100000 


0.255569 


0.21104 


1.9549 


102.65 


451.90 


0.000316228 


0.255847 


0.21156 


1.9539 


102.23 


448.34 


0.00100000 


0.256701 


0.21321 


1.9518 


100.94 


436.59 


0.00316228 


0.259029 


0.21814 


1.9488 


97.304 


399.80 


0.0100000 


0.267020 


0.23387 


1.9269 


86.717 


317.83 


0.0316228 


0.287014 


0.27830 


1.8859 


65.176 


173.48 


0.100000 


0.331664 


0.39429 


1.8758 


35.484 


49.785 


0.316228 


0.428943 


0.71809 


1.8628 


13.582 


3.6035 


1.00000 


0.771803 


1.6210 


2.2704 


2.8947 


0.0076 


3.16228 


0.984801 


4.7134 


2.4773 


0.34516 


0.0000 


10.0000 


1.26682 


14.999 


2.4199 


0.03598 


0.0000 


31.6228 


1.60590 


48.113 


2.3283 


0.00371 


0.0000 



Table 3 

The results for the regularized quark condensate and for the regularized quark dis- 
persion relation obtained with the ansatz of Eq. (46) . All results are in dimensionless 
units ofcj = 0.19 GeV^ = 1. 



of the hadron spectra, it is convenient to write the dispersion relation as a 
sum of the analytical infrared term of Eq. (44), plus the free quark dispersion 
relation dominating the ultraviolet, and plus a finite and compact term E{p), 
an integral that we compute numerically. 



E[p) = 7, 7, T-r^ 

Ijj, n p"^ — m[pY 



+ pC{p) + moS{p) + E{p) 



(45) 



The numerical integral E{p) decays like 1/p^ in the chiral limit of small current 
quark masses and decays like 1/p^ in the case of large current quark masses. 
E{p) is negative and we can conveniently fit it with the rational function, or 
Pade approximant with odd powers of p only. 



E{p) 



I 



eo + eip + e-ip-^ + e^p^ 



(46) 



We show the best fitting parameters Cq, ei, 63, 65 in Table 3. With our excellent 
fits of the dynamical quark mass m{p) and of the quark dispersion relation 
E{p) we are well equipped to address further problems, such as the breaking 
of chiral symmetry or the hadronic excited spectra at finite temperature T. 
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Fig. 5. (a) We show a log log plot of minus the regularized quark condensate 
— {tpt/j) + {ipilj)o as a function of the current quark mass mo, with vertical dot- 
dashed lines representing the current masses of the quarks u, d, s, c, b, t. (b) We 
represent the regularized quark dispersion relation E{p) — ^ with increasing num- 
ber of dashes per curve for five different current quark masses tuq with respective 
values 10""^, 10~^, 10~^'^, lO'', 10^'^. The quark condensate has an infiection point 
for finite quark masses close to the strange quark mass, and for large masses it rises 
linearly with uiq. All results are in dimensionless units of cr = 0.19 GeV^ = 1. 

5 Conclusion 



While the chiral limit of rrio << ^/a was already well known in the literature, 
we find unanticipated effects for finite mo — \/o" and for heavy nio » y/a 
current quark masses. We study in detail the large tuq limit, performing an one- 
loop expansion in the dimensionless number a/niQ^. We also develop a new 
technical approach to solve the mass gap equation, utilizing the variational 
principle to increase the precision of our mass solution m{p) in the infrared 
limit of p — )■ 0, relevant to compute the mass gap m(0), an order parameter 
for the chiral phase transition. We also show that the dynamical generated 
constituent quark mass mij)) can be quite well fitted by our inverse even 
quartic polynomial ansatz, a Pade approximant with parameters co, C2 and 
C4 depending only on the current mass niQ. 

Our surprising results are that the dynamical mass generation has finite effects 
persistent beyond the chiral limit. At rriQ ^ ^/a , in particular for masses 
similar to the strange quark s mass, the quark mass generation m(0) — rriQ 
is maximum, as shown in Fig. 4, while one would naively expect the quark 
mass generation to be maximum in the chiral limit close to the up u or down 
d quark masses. A second order parameter, the regularized quark condensate 
{ipij) ~ {'^'^)o, monotonously grows in absolute value with rriQ and shows an 
inflexion point for masses similar to the strange quark s mass, as depicted in 
Fig. 5. In the limit of heavy quark masses rriQ » ^Ja of the charm c, bottom 
h and top t quarks, although the mass gap difference m(0) — mg vanishes like 
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cr/?7io, it occurs that the energy density difference £^ — £^o is maximum. In the 
heavy quark hmit, the energy density difference converges to a constant hmit 
when rriQ — )■ oo, and we show in Fig. 2 that it is one order of magnitude larger 
than in the chiral hmit. This may be relevant for cosmology, contributing to 
the dark energy. 

The numerical variational technique developed here, together with the detailed 
solutions for the running quark mass m(p) and for the quark dispersion relation 
E{p) as a function of the current quark mass itlq and of the string tension a, 
are necessary tools for the continuation of our program to study the QCD 
phase diagram and the hadron spectrum at finite temperature T and density 
/i, when the string tension a becomes quite small, possibly smaller than the 
light current quark masses. 
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